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Abstract 

Data on the total cross sections for the pp PPVi J"^ ~^ pnr], and 
pn — > drj reactions and the pp — > ppr] differential cross section near thresh- 
old are analysed in a one-meson-exchange model. After including initial 
and final-state nucleon-nucleon distortion, the magnitude and most of the 
energy dependence are well reproduced. It is found that the contribu- 
tion of /9-exchange is larger than that of 7r-exchange. With destructive 
p/tt interference in the pp case, the model explains quantitatively the 
pp pprj/pn pnr] cross section ratio and the slope of the pp — > ppij 
differential cross section. Such an agreement would be destroyed by any 
significant r/-exchange term. The residual energy dependence may be as- 
sociated with ry-nucleon rescattering that has not been taken into account. 
The pn — > pntj/pn dt] ratio depends weakly upon the nature of the 
particle exchanges, being determined primarily by the nucleon-nucleon fi- 
nal state interactions. The proton analysing power is predicted to remain 
small in the low energy region. 
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1 Introduction 



The experimental database on rj production in nucleon-nucleon scattering has 
expanded significantly in recent years. The excitation function of the pp —>■ 
ppr] total cross section shows such a rapid rise with excess energy Q ||, 
^, that the major error is often associated with the determination of Q rather 
than of the cross section itself. Equally striking is the large cross section found 
for quasi- free t] production on a deuterium target under conditions of well- 
controlled kinematics the pn pnrj/pp ppr] ratio is found to be over 6 near 
threshold ^ . The CELSIUS group has also measured the pn dr] total cross 
section at well-defined cm. energies ||^, 0] and found the first evidence of A^*(1535) 
dominance from the energy dependence. More tantalisingly, they showed that at 
very small Q there appears to be a threshold enhancement which, though not 
as spectacular as that found in the analysis of early Saclay data, nevertheless 
indicates a strong ?7-deuteron scattering length. 

The only measurement of the angular distribution of r^-mesons near thresh- 
old suggests that (i- waves cannot be neglected in the ri-{NN) system for Q ^ 



30 MeV [IC]. It is here important to note that the differential cross section has 
the opposite curvature to that of n^p —>■ rjn and we will argue that this helps 
to identify the reaction mechanism. In view of these new data, as well as the 
experimental advances in tt"p rjn and especially 'yp —>■ rip, we believe that it 
is useful to revisit the phenomenological analysis of these reactions. 

In the usual theoretical approach, the A^*(1535) (or other) isobar is excited in 
nucleon-nucleon collisions through the exchange of a single meson X |TT], 0, 113 



14, |T^, 16, 17|. The A^*(1535) has a large branching ratio into t] N so that, after 
its decay, one is left with an t] plus two nucleons in the final state. Though the 
philosophies are generally rather similar, the theoretical calculations differ in their 
details as to which exchanges are relevant, the structure of their coupling, and the 
importance of the associated form factors. Different techniques have been used in 
order to include the effects of the initial nucleon-nucleon distortion and the final- 
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state interactions. The neutron-proton final-state interaction can also produce a 
deuteron in the exit channel and, as we shall see, the ratio of the pn drj and 
pn npr] cross sections is determined primarily by low-energy neutron-proton 
dynamics with a weaker dependence upon the production mechanism 

It is generally agreed that in the near-threshold region, Q < 40 MeV, the 
energy variation of the total cross section for the NNrj channel is fixed mainly by 
the factor coming from phase space, modified by the nucleon-nucleon final- 
state interaction. It is known, however, that rj rescattering does produce threshold 
enhancements in the pd ^Her] and dd — ^ '^B.er] reactions and it is 
likely that the residual energy dependence in the pp pprj excitation function 
could be due to such an effect. The energy dependence coming from the input 
single- nucleon production amplitudes is rather modest over a small range in Q. 

The relativistic Born amplitudes for single-meson exchange are first evaluated 
at threshold by neglecting all distortions. To include the effects of the final-state 
interactions it is, however, necessary to construct a potential from the amplitudes 
and this is most easily done in configuration space. The energy and angular 
dependence coming from vr^p —* rjn and the other elementary amplitudes can 
then be introduced in a perturbation approach. 

The kinematics of the processes are outlined in Section 2. As exchanged 
particles, we consider the vr, rj, p, and u. The forms of their coupling to the 
nucleon are given in Section 3, together with the coupling constants used in this 
work, though it must be stressed that the uncertainty in the rjNN value is very 
large. Values for the XN —>■ rj N amplitudes are also discussed in this section. 
There are, of course, measurements of the pion-induced cross section and the 
r^-nucleon elastic cross section can be deduced from this inside a unitary model. 
To obtain the necessary vector meson information, we interpret photoproduction 
data within the framework of the vector meson dominance approach. 

The formulae for the bare meson exchange A^A^ NNrj and pn drj ampli- 
tudes are derived in Section 4; the energy and angular dependence coming from 
the XN T] N amplitudes are there introduced. 
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Section 5 is devoted to the evaluation of the final state interactions in the pp 
and np systems. It is found that the enhancement factor varies very fast with 
excitation energy, and also depends upon isospin in the A^A^ channel as well as 
on the mass of the exchanged meson. In contrast, the distortion of the initial 
nucleon-nucleon wave is slowly varying and we try to include its effects simply 
by using the imaginary part of the relevant phase shift. The //-nucleon final- 
state interaction is not taken into account, since it would require a consistent 
three-body treatment in order to include this simultaneously with the nucleon- 
nucleon interaction. There is, moreover, as yet no credible ?7-nucleon potential 
from which to calculate the wave function at short distances. The two-body 
pn dr] amplitude is discussed in Section 6. 

As shown in Section 7, p-exchange is more important than 7r-exchange, with 
the uj giving only minor modifications. Taking the p/vr interference to be de- 
structive in the proton-proton case, these three exchanges reproduce well the 
measured 6.5 : 1 pn : pp ratio. Since pure 7] exchange would lead to a ratio of 
below 1, the results yield an upper bound on the poorly determined r^-nucleon 
coupling constant. With standard parameters, the model reproduces well the 
magnitude and most of the energy dependence observed in the total cross. There 
is an indication that the model underpredicts the pp pprj data at low Q and 
the same is true for the two-body pn drj results. This is probably due to the 
neglect of rj rescattering. On the other hand, the reasonable agreement found for 
the a{pn — > pnvj) / a{pn drf) cross section ratio is mainly a refiection of the np 
final state interaction. 

The dominance of the p and the destructive effect of the vr explains also the 
shape of the rj angular distribution in pp —>■ pprj, that is larger at 90° than in 
the forward direction \TU\. There are still insufficient input data to give definitive 
predictions for the proton analysing power, but the indications are that this 
should be rather small in the near-threshold region. Our conclusions are given in 
Section 8. 
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2 Kinematics 



We wish to describe the pp — > pprj, pn — > nprj and pn — > drj reactions in the 
near-threshold domain. The energies of the incident nucleons in the cm. system 
are denoted by E and their momenta ±p. At threshold they are related to the 
nucleon and rj masses, M and /i, by 

E = M+^n 

p = V/^M+i/i2. (2.1) 

The excess energy Q is fixed by the total cm. energy as 1^ = 2E — 2M + 
rrir, + Q, where we shall neglect the neutron-proton mass difference except in the 
determination of Q from the beam energy. In the case of the two-body pn drj 

reaction, the 2M is replaced by the deuteron mass, M^. Reduced masses in the 
final state are defined by 

As an alternative to Q, experimentalists often quote the value of the maximum 
Tj momentum, p^, in units of the rj mass; 



/i V 

If the matrix element A4 were constant at its threshold value, the only energy 
dependence in the unpolarised total cross section would come from the factor 
arising from phase space. Thus, for the reaction pp — > ppr), 

where the summation is over the initial and final spin projections. This expres- 
sion incorporates a factor of one-half coming from the identity of the final-state 
protons. Such a factor is absent from the corresponding pn — > pnr) formula. 
Finally, we also need the cross section for pn ^ drj, which is 
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As we shall see later, the strong nucleon-nucleon final state interaction will 
lead to a major modification in the behaviour such that it will no longer be 
permissible to factorise the dynamics from the phase space, as has been done in 
Eq. (0). 
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3 Input Amplitudes and Vertices 



The basic inputs required for the evaluation of our Feynman diagrams are vertex 
functions and r^-production amphtudes Ai, which are the matrix elements of 
iCint = —iTi-int- We list below their general forms as well as the corresponding 
approximations used in the evaluation of the near-threshold amplitudes. 

3.1 Pseudoscalar meson-nucleon vertex 

The standard nNN coupling in terms of four-component spinors u{p) is 

M^NN = -— uip'Hif- ^) r ■ ct>Mp) , (3.6) 
where we take the pion-nucleon coupling 

^ '-^ (3.7) 
to have the numerical value G'^/Att = 13.6 Since we are specialising to near- 



threshold reactions, it is sufficient in our applications to put p' = and E' = M 
in the final state. 

It is important to take into account the off-shell extrapolation in the vertex 
function and for this we introduce a monopole form factor 

FAq') = ^^f^ , (3.8) 

where is the square of the pion four-momentum. We take the value A^^ = 
1.72 GeV/c, as recommended in pl| . 



The size of the ?7-nucleon coupling constant is extremely uncertain, with values 
of G^/4:TT between and 7 being quoted in the literature (see e.g. [0). A value 
of 0.4 has been deduced from fits to photoproduction [^. However, even with 
a value as small as that, it is claimed that T^-exchange would strongly infiuence 
7] production in nucleon-nucleon collisions [T^. Since pure r^-exchange gives a 
completely wrong prediction for the ratio of the relative production in the pp and 
pn cases, these data will provide evidence for some upper limit for G^ / Air within 
the framework of the meson-exchange model for rj production. 



The range parameter in the rj form factor will be taken, as in the Bonn meson- 
exchange potential to have the value = 1.5 GeV. 



3.2 Vector-meson- nucleon vertex 

The generic form of the couplings of vector mesons to the nucleon is (see e.g. 



MvNN = -igv u{p'. 



= -igv^-V, (3.9) 

where k = p' — p is the four-momentum of the vector meson and e = e{k) its 
polarisation four-vector. 



The numerical values of the p and u coupling constants are, according to 

gp = 3.25, Kp = 6.1, , . 

^^ = 15.9, K^ = 0.0. 



As in the Bonn potential ||2^, ^ we shall use the monopole form factors of 



Eq. (|3.8| ) with numerical values Ap = 1.4 GeV/c and A^^ = 1.5 GeV/c. 

To take account of the 1 = 1 nature of the p meson, the vertex function 
in Eq. ( p.9[ ) must be multiplied by an isospin factor t ■ <pp in analogy to the 
pion-nucleon vertex. 

3.3 Pseudoscalar meson-induced rj production 

In the absence of polarisation data, the spin-non-flip vrA^ rjN amplitude is 
sufficient to describe the bulk of the experimental data near threshold. In the 
two-component reduction this has the structure 

MinN ^r]N) = u{p')[-iK] t ■ (j)^ u{p) , (3.11) 

where /i^ may depend upon the production angle 6. For neutral pions the cm. 
differential cross section becomes 

whereas for charged pions the cross section is a factor of two bigger. 



Until the large Crystal Ball data set is fully analysed , the available results 



on the IT p ^ rjn differential cross section are generally far less systematic than 
those for photoproduction. Parameterising the differential cross section as 

^ = j [feoi^olcos^) + 61 Pi (cos ^) + 62 P2(cos0)] , (3.13) 



the threshold value was determined many years ago as Awba = (9.1 ±0.8) mb [ p7[] . 
This is a little higher than the first Crystal Ball measurement at 720 MeV/c, that 



gives AttBo = (7.4 ±0.3) mb |2^. This group has also a preliminary determination 
of the shape of the differential cross section at (716 ± 16) MeV/c, which leads to 
bi/bo = —0.042 and &2/&0 = 0.137 [^. Assuming an energy dependence of the 
total cross section similar that for photoproduction shown in Eq. ( |3.26|) , we fix 
the n^p — > rjp amplitude squared to be 

IKl"^ = [(128 - 3O3r/2)Po(cos0) ± 23r7Pi(cos^) ± 3427^2 P2 (cos ^)] mb/sr . (3.14) 

There are as yet no measurements of the proton analysing power near thresh- 
old, but the smallness of the Pi coefficient in Eq. (|3.14| ) suggests the angular 
dependence arises primarily from the interference of a small spin-non-flip rf-wave 
with the dominant s-wave. This could for example be due to the tail of the D13 
resonance. We therefore neglect the Pi term and take 

K = K + \hl (3(p, ■ kf - 1) , (3.15) 



where 

|2 = (128 - 3037]^) mb/sr , Re[hl/hl] = 1.3ri'^ . (3.16) 

Multiresonance fits to projections of 7r~p —>■ rjn and ir^p TT~p/TT^n data 
onto the 6*11 channel allow one to estimate the //-nucleon elastic scattering am- 
plitude. Although the results, especially on the real part, are rather sensitive 
to the assumptions made |]30|, 0, a typical value of the scattering length is 
(^r^N = (0.83 ± 0.27i) fm. The elastic ?]-nucleon scattering amplitude correspond- 
ing to Eq. (|3TTD is 

h^ = -An(^l + -^^ a,iv. (3.17) 
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The exchanged mesons are far from the mass shell in our model. We assume 
that the effects of this can be taken into account by multiplying the pion-exchange 
amplitude by the monopole form factor FT^{q^) of Eq. (3^), and similarly for rj 
exchange. 

3.4 Vector-meson-induced 77 production 

Although nothing is known directly about 77 production by vector mesons, one 
can use the vector-meson dominance model (VMD) to relate photon- and vector- 
meson-induced reactions [p^ . 

The isovector component of the photon couples to the p field and the isoscalar 
part to the tu, though such an identification is only valid for transversely polarised 
vector mesons. Neglecting higher mass vector mesons, the photoproduction am- 
plitude may be written in terms of off-shell strong-interaction amplitudes as 

6 6 

M{'^p ^ r]p) = — M±{poP -> VP) + -^M±{ujp r]p) , (3.18) 

with an analogous relation for neutron targets. 

Since universality is broken in the VMD model, we choose as effective p-meson 
coupling constant the geometric mean of the two values derived by Benayoun et 
al. H, 

^ = 2.4. (3.19) 
47r 

This number is in reasonable agreement with those extracted from p-photoproduction pSj . 

The value of the uj coupling constant is less certain. Here we have chosen to 
scale it to the p coupling according to the r(p — > e+e~) and Vioj e+e^) decay 
widths, giving 

^ = 27 . (3.20) 

Data on the photoproduction of ?7's on deuterium show that near threshold 
(T(7n — i> rin)/a{'-)p rjp) = 0.66 ± 0.03, where some attempt has been made to 



include the systematic error arising from the use a deuteron target |32|. The low 



cross section found for the coherent "yd — rjd reaction implies that isovector 
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photons dominate the production, so that the amphtude ratio is negative: 

r = — — = 0.81 ± 0.02 . 

jvl ['-fp —>■ r]pj 

The vector-meson-induced amphtudes thus become 



(3.21) 



/ /I + ^ 



J^±{ujp — i> rip) 



e V 2 
fu; fl-r 

e 



M (7P rjp) , 
M. {-yp Tjp) . 



(3.22) 



The threshold photoproduction amphtude is proportional to cr-e^. The gauge- 
invariant photoproduction amplitude that reduces to the correct threshold S-wave 
limit is 



Mj = R^ |m(p')75(/^- f i)u[p) — ^Svrl^Eo+T^V-eC 



where 



{E + M){E' + M){W - M)R 



7 ; 



(3.23) 



(3.24) 



with k the photon momentum four-vector and W the cm. energy. 

Assuming, as for the pion-induced reaction, that there is also a contribution 
from the D13 resonance as well as the Su at low energies, the more general 
reduction is 



W T]UEQ+rT ■ e + \E2-\'i(T ■ q q ■ e - (T ■ e\) ( 



(3.25) 



Taking the (i-wave to first order, the unpolarised photoproduction differential 
cross section is related to these multipoles through 



^ilP - VP) = y l^o+r (1 + ReiE^jEo 



Siq-kf 



(3.26) 



The near-threshold 7p rjp differential cross section can be parame 



terised ||3^ 
da k 



dVL k 



[(4.59 - 10.9 7]^) Po(cos^) - O.291r/Pi(cos0) - 3.2 V P2(cos 0)] /ib/sr . 

(3.27) 
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Neglecting the Pi term, which must arise from an interference with a small p-wave 
component, we deduce that 



l-£'0^ 



(4.59 - 10.9 /xb/sr , i?e(£'2-/S, 



-0.70r/^ . 



(3.28) 



At threshold \Eo+\ = 0.0214 fm and hence \R^\ = 0.146 fm^. The energy variation 
of Eq^ is mainly a reflection of the dominance of the A^*(1535) resonance near 
the 7] threshold. 



In the same approximation, the proton analysing power is p6 
Ay = —3 Jm(i?2-/-E'o+) sin 9 cos 9 . 



(3.29) 



Ay remains small for photon energies below 850 MeV |3^, with evidence for both 
sp and sd interference. This gives 



Im{E2-/Eo 



-3 ± 0.5)r/2 + (16 ± 4)r/^ 



(3.30) 



We now make the ad hoc assumption that all three polarisation states have 
the same strength in the vector-meson-induced reaction. The amplitude is then 
proportional to cr-ey, with the proportionality constant being determined by the 
vector-meson-dominance model. The relativistic form which has this limit is 



M{Vp r]p) = Rvu{p') Ih^i/ 1^-1^/)+ mv'j5 



u[p) 



{E + M){E' + M){W-M + mv)Rvri^o- • eC .(3.31) 



The scale factors Ry are determined from the VMD model, Eqs. ( p.22| ), evaluated 
at the 7] threshold W = M + fi: 



p _ p /p 

e 



fx 



Ru) — R-. 



(3.32) 



From the value of R^, we deduce that \Rp 



e \ 2 

0.916 fm^ and \RJ = 0.310 fm^ at 
threshold. 

As before, we attempt to describe the off-shell extrapolation through a form 
factor. Since we are starting from data with zero-mass photons, the correct form 
factor at the p-induced photoproduction vertex is 

W) = Fp{q^)/FM ■ (3.33) 
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We take for Fp{q'^) the monopole form factor of Eq. ( ^.81 ). This is not in contra- 
diction with the off-shell behaviour of the ?7-photoproduction amplitude measured 
via electron scattering at modest |^ . 



3.5 Deuteron vertex 

Both nucleons are close to their mass shells at the deuteron vertex, so that each 
pole gives a contribution after integration over the internal energy variable. The 
deuteron vertex matrix element corresponding to the deuteron S- and D-state 
wave functions 05 (r) and (poi^) is 
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X 



iSUk){-ir-e*)iSF{pa-k)^ 
■ 4 MQr) + {o- • 4 - • Oil 4 • Qr} (I>d{Qr)] C,u{Pd-k) 



2M 



{27r){S{k'-En)+S{p',-k''-E,)}. 



(3.34) 



En and Ep are the nucleoli on-shell energies so that, e.g. E^ 
Qji is the proton-neutron relative momentum. 



k'^ + M^, and 
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4 Amplitude for pp — > pprj 



We describe the reaction pp pprj in terms of Feynman diagrams, the prototype 
of which is displayed in Fig. 1, considering the four exchanges, p, uj, tc, and r/. 
The matrix elements are evaluated at threshold, where the momenta in the final 
state vanish and the values of the corresponding incident proton energies and 
momenta are given in Eq. ( |2.1| ). We first study neutral-particle exchanges; the 
necessary isospin factors will be derived at the end of the section. 



Pn 



V 



P 



Q-2Pv 



a) 



7r,p 



Q-2Pv 



T 



k 



k 



-i- 



^^^V I 1 

^ -1- 2Pr, 



P 



q - 



-P 



Pv 



Figure 1: Meson exchange diagrams for pp —>■ pprj. Diagrams c) and d) are 
obtained from a) and b) by interchanging the final state nucleons. 



4.1 Vector-exchange diagrams 

The propagator for vector meson exchange is 



tD,Ak,m) = '^^f±^ . (4.35) 



m 



V 



When the vector meson four-momentum fc^ is contracted into the VNN vertex, 
it vanishes identically. We thus obtain for exchange in diagram a 

— 1 

Ma{pp PPV) = M^{p pp°)7^ :rM^{p°p rip) . (4.36) 



k'^ — m 



Using the expressions for the matrix elements of the vertex ppp and amplitude 
pp rjp, given in Eqs. ( |3.9|) and ( p.30| ), leads to the threshold expression 

Maipp - PPv) = f"^" ^ [f2(-Mp)/F,(0)1 JCa , (4.37) 

Mp + nip L ^ -I 
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with 

Ka = a [7^1 <T ■ p Ci] [4 C2] - ih [4 o- Ci] ■ [4 P X a C2] . (4.38) 
The parameters a and h are given by 

a = (M-m,)(l-^), 

h = |(mp + ^)(l + /tp) . (4.39) 

The momentum of the intermediate vector meson is = —p and the form factor 
is evaluated at mass squared k"^ = —Mfi. 

We can similarly write down the expression for diagram b, where the two 
vertices are interchanged: 

JCb = -a [v^ Ci] [vl (T■pC2]+^b [4 P x o" Ci] ■ [vl o- C2] , (4.40) 
with k = +p. 

It is easy to interpret the 6-terms if we explicitly write the sum over the 
polarisation states of the p meson. Thus, in Eq. ([4.38|) , this gives 

[V3 o- Ci] [^4 P X o- C2] = E IVs ^ ■ ep Ci] [^4 ep ■ (P X cr) C2] (4.41) 

Here we recognise the vertex for the pNN coupling as ■ (k x cr) and that for 
p-induced 1] production as cr ■ . The small a-term arises from the part of the 
pNN coupling that contains the vertex factor Cp ■ k . 

At threshold the initial total spin of the pp system is = 1 with = in the 
final state and we shall let e denote the spin vector of the spin-one pp state and 
7] the spin-zero state vector. The contribution of diagrams a and b is the same as 
that for the crossed diagrams c and d: 

ICa + ICb = }C, + }Cd = -2(26 -a)p-e,7]f. (4.42) 

From Eq. ( [4.37] ) and ([4.42|) , the total p° contribution becomes 



Mp{pp PPT]) = {Ap - 2Bp) p ■ Si r]f , (4.43) 
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where 



•^0 



16M 



16M 



Mfi + ml 



9pRp 



f;(-m/.)/f,(o) 

F^(-M^)/F,(0) 



b. 



(4.44) 



With the values of the parameters given in section 3, it can be seen that the 
^-term is of no numerical importance for p exchange. 

When the ci-waves are included in the dominant 6-term, the amplitude of 
Eq. ( ^.43| ) is modified slightly to read 

Mp{pp pprj) = {Ap-2Bp) p ■ Ei + \{E2- / Eq+){?>p ■ q q ■ Ei - p ■ Ei) rjf . 

(4.45) 

To first order in E2-/E0+, this gives exactly the same form as in the 'jp rjp 
cross section of Eq. ( |3.26| ). 

The amphtudes for u are completely analogous to those for p exchange, so 
that at threshold 



M^{pp PPT]) = (X - 2B^) P- Ei7]f 



(4.46) 



where the amplitudes are as defined in Eq. ( [4.44|) and the parameters in Eq. (|4.39|) . 



4.2 Pseudoscalar-exchange diagrams 

The four diagrams with 7r° exchange are labelled as in Fig. 1. The contribution 
from diagram a is 



Ma{pp — * ppv) = M{p pn^) 



^2 _ ^2 



■ Miir^p rip) . (4.47) 



Inserting the forms for the npp vertex and r^-production amplitude, given by 
Eq. {^) and (gJ^, into Eq. (CT) yields 

1 



Maipp-^ PPV) = 2MG^ K 



where 



Mp + 



Fj{-Mp) Ka 



(4.48) 



(4.49) 
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Using the same techniques as for p exchange, the sum of the amphtudes for 
diagrams c and d is the same as that for a and h: 

]Ca + }Cb = }Cc + ICd = 2p-e,r]f. (4.50) 

The final result is then 

M^ipp ppr]) =1)^ p- Ei r]f , (4.51) 

where 

= 8MG^ K 



2 L^.'(-M/i) . (4.52) 

With the inclusion of rf-waves, the input amphtude /i^ has an angular dependence 
given by Eq. ( |3.15| ). 

The expressions for the r^-exchange amplitudes are algebraically identical to 
those of pion exchange, with the vr index being replaced by an t]. Of course, as 
discussed in Section 2, poor knowledge of the value of the t]NN coupling constant 
makes the strength of any rj contribution extremely uncertain. 

4.3 Isospin considerations 

The expressions derived for the various amplitudes correspond to neutral-meson 
exchange for the pp case. The NN NNt] isospin-one production amplitude is 

Ml = VfP- [{A, - 2Bp) + (X - 2B^) +V^ + P,] x) ■ X^ ■ (4-53) 



It is straightforward to generalise this to the case of isospin zero [11 



Mo=p- e} 7], [-3{Ap + 2Bp) + {A^ + 2B^) - 3V^ + V^] <p] <\>, . (4.54) 

Here x and are isospin-1 and -0 operators and the subscripts refer to the 
particles exchanged. 

The spin-average of the square of the matrix element becomes 

1 ^ \M{pp ^ pp^t = ^^\{A,~ 2B,) + {A^ - 2B^) +V^ + V,\' (4.55) 
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for the pp case and 

2 

lY.\Mpn^nprj)f = ^ [|-3(^, + + (X + - + 

+ \{Ap-2Bp) + {A.-2B^)+V^ + Vr,f] (4.56) 
for pn. Application of the formulae in Section 2 then give the cross sections. 
4.4 Kinematics 

Since we are only considering cases where the excess energy is low, the transfor- 
mation between the final rj N and rj NN systems is non-relativistic. Taking the 
laboratory momentum to be the same in the 77 N and 77 NN cases, one sees that 
the cm. momentum is lower in the single-nucleon than the two-nucleon case: 

fl + a/2M\ 
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5 Initial and Final Nucleon-Nucleon Distortion 



Most of the energy dependence observed in the cross sections for the A^A^ — *■ 
NNf] reactions near threshold can be ascribed to the behaviour of the three- 
body phase space that has been modified by the very strong nucleon-nucleon 
final state interaction (fsi). It has been traditional to incorporate the effects of 
such an fsi by multiplying the predicted cross section by a Watson enhancement 
factor. This prescription gives a rapid energy variation parameterised in terms 



of the A^A^ scattering length and effective range It does, however, give too 
steep a fall at higher A^A^ relative momenta 0] and does not attempt to provide 
an overall normalisation factor. 

In their analysis of the experimental data, the authors of [|ry] took the en- 
hancement factor as the ratio of the squares of the interacting A^A^ S-state wave 
function to the equivalent plane wave. In the absence of a model for the re- 
production operator, they evaluated the ratio at a fixed radius of r = 1 fm, 
which is close to the maximum of the A^A^ density |^ . 



To obtain more realistic enhancement factors, we must study further the 
meson-production operator. In the plane-wave approximation, the dominant p- 
exchange contribution from the diagram of Fig. la can be described by a potential 

i + -rp^- (5-58) 

m*r J m*pr 

In the vicinity of the rj threshold, the energy transfer is shared equally by the 
two nucleons, so that the range of this propagation is determined by the reduced 
p-mass, given by 

mf = ml- uj^/A , (5.59) 

where uj is the total energy of the 1] meson. 

Taking matrix elements of the potential between interacting nucleon but t] 
plane waves and projecting out the threshold angular momentum transition, we 
are left with the amplitude 

M';' = MpCp ^V\r) ri(m;, r) jo(br,r) V-^+^r) dr , (5.60) 
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with M.p the appropriate plane wave p-exchange amphtude of Eq. ( |4.53| ) or ( |4.54| ), 
and where 

Cp = ^m;\m;^+p^) (5.61) 

and 

("1 \ — Tfl* V 

1 + ^^- (5-62) 
m*r J rripr 

The wave function ipkif ) describes the S'-wave of the final / = !, J = Oor/ = 0, 
J = 1 NN system and ipp{r) the incident high energy A^A^ P- waves. 

The exchange of other negative parity mesons, {u, vr, 77), leads to a very similar 
structure, though with different effective masses given by Eq. ( |5.59|) . 



When monopole form factors of the type given in Eq. ( |3.8| ) are introduced 
at the pNN and production vertices, the only modification to the formula is to 
replace the Fi(m*,r) propagator in Eq. (|5.62|) by 

* , (A*2+p2)2 1 1 /2 2A 



(A* 2 - m*2)2 
with 



(5.63) 



A*2 = a2-c^74 (5.64) 

For r ^ 1/A*, Y\{m*^^r^ = 0{r) and so the form factor reduces a little the 
sensitivity to the small-r components of the wave function. 

In order to investigate purely the effect of the A^A^ fsi near threshold, ipki^) is 
replaced by the plane wave spherical Bessel function ji{pr). Near threshold the 
7] plane wave factor jo{^Prir) can be taken to be unity. We then quantify the final 
state enhancement through the ratio of the amplitudes 

E,{k) = M'^'/Mr (5.65) 

calculated from Eq. (|5.60|) with interacting and plane A^A^ final-state waves. 

The enhancement factors l-E^P are shown in Fig. 2 for both / = and 1=1 
NN states in the case of p-exchange. It can be seen that Coulomb effects are 
significant for pp final states with k < 0.2 fm~^. For vr-exchange the energy de- 
pendence is almost identical, though the overall magnitude is somewhat less. To 
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illustrate this, the 1 = enhancement factor for vr-exchange is shown multiplied 
by a factor of 1.85. 

Because the pole of the antibound state in the pp system is closer to threshold 
than that of the deuteron in np, there is significant energy variation in the ratio 
of the {1=1)/ {1=0) enhancement factors. Nevertheless, as can be seen in Fig. 3, 
for fc^ > 0.2 fm~^ the ratio is roughly constant at 1.85. This deviation from 
unity is important in the understanding of the pp — pprj/pn npt] ratio. The 



empirical approach used in |T0|, where the wave function squared is evaluated at 
a fixed radius of 1 fm, would lead to a ratio with a similar energy dependence to 
that shown in Fig. 3 but lower in magnitude by an overall factor of 0.8 




Figure 2: Enhancement factors, defined by Eq. (|5.65|) , evaluated for isospin-one 
and -zero NN S"- waves using Paris wave functions |^ . The solid lines are without 
the Coulomb repulsion, whereas the broken line includes this effect in the 7 = 1 
case. The calculations were done for p-exchange but vr-exchange leads to a very 
similar energy dependence though reduced in magnitude by a factor of about 1.85 
in the 7 = case and 1.65 for 7 = 1. This is illustrated for 7 = by the dotted 
curve that has been scaled up by 1.85. 
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Though the initial-state interaction is expected to vary httle with energy, 
the evaluation of an A^A^ wave function from a potential at energies as high as 
1300-1400 MeV is extremely dubious. In order to take the initial flux damping 
into account, the amplitudes are multiplied by the factor r]^ = e~^''^^^'^\ This is 
typically 0.77 for the ^Pq state and 0.73 for ^Pi. 




1 I I I I I I 

0.2 0.4 0.6 0.8 1 



Figure 3: Ratio of enhancement factors for isospin-one and -zero NN final states 
evaluated for p-exchange using Paris wave functions with Coulomb effects 



The ratio of the squares of the wave functions at r 1 fm is similar in shape but 
reduced in magnitude by about 20%. 



We have neglected the distortion of the final rj wave. This is, in part, due to the 
difficulty of including consistently final state interactions simultaneously in the pp 
and rip channels. There is, furthermore, a lack of a realistic ?7-nucleon potential 
in the literature. On the basis of potential-model fits to ?7-nucleon scattering 
parameters, Garcilazo and Pena predict a strong threshold enhancement in 77- 



deuteron elastic scattering However, probably as an artefact of their one- 

term separable assumption, their potentials are extremely attractive at short 
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distances (r < 0.2 fm). Since 77-production is sensitive to wave functions at small 
r, using these potentials as a basis for estimating final state interaction effects can 
lead to very unrealistic values. We have therefore represented the rj hj a plane 
wave though, as we shall see later, there are indications from the production data 
that this is inadequate. 
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6 Amplitude for pn — > drj 



The triangle diagrams describing the pn — > drj reaction are shown in Fig. 4. They 
are very similar to those for pn — > nprj with the addition of a neutron-proton 
final-state interaction to produce the deuteron. 



y Pr, 

Pi y Pi Pi Pd 




\ Pr, 



Figure 4: Meson exchange diagrams for pn ^ drj. 



The Feynman amplitude corresponding to Fig. 4a for p exchange is 

X tS'^M) (-^r ■ e^) zSf{p[) M>'{p'>p ^ VP) u{pi) , (6.66) 

where we have treated the incident neutron in the charge conjugate representa- 
tion. Here k is the four-momentum of the p-meson and we can neglect the term 
proportional to k^ki, in the p propagator since this gives only binding energy 
contributions. The isospin factor of 3/^/2 takes into account the contribution of 
charged intermediate mesons. There is, in addition, the diagram of Fig. 4b where 
the 7] production takes place at the neutron vertex. 

We calculate the triangle diagram in the spectator approximation by perform- 
ing the integration over after suitably closing the contour of integration. This 
means exploiting formula of Eq. ( ^.34| ). After some straightforward manipula- 
tions we get in configuration space, as the sum of both p-exchange diagrams, 
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X Cp [i^sir) + 2^z)(r)] Fi(m;, r) jo(|p„r) V^W(r) rfr . (6.67) 

where yi(m*,r) and Cp are defined in Eqs. ( p.63| ) and (|5.61|) . 

In a similar fashion, we find the pion-exchange amplitude to be 

X a / [^5(r) + 2V^B(r)] Fi(m;,r)jo(b.r)^(+)(r) (6.68) 

with Fi(m*,r) and being defined as in Eqs. (|5.63|) and (^.61|) , but with nip 
replaced by m^. 

The contributions from uo and r] exchange can be derived immediately from 
Eqs. ( |6.67|) and (|6.68| ) in an analogous manner to that for pp pprj. 

The ?7-production operator, especially that part corresponding to p exchange, 
is of short range. At short distances and low energies the deuteron S-state wave 
function and the / = scattering wave functions at low energies are related by 



the extrapolation theorem 



^'C\r) ^ -^=^=— i^sir) e''^ , (6.69) 



where 6s is the 5- wave triplet phase shift. If we neglect the S-D coupling, the 
relation becomes exact as fc^ —o? = —M Ed, where Ed is the deuteron binding 
energy. 

By comparing Eq. ( |6.68| ) with Eq. (|5.61| ) , corrected for isospin and form factor 



effects, the extrapolation theorem allows us to estimate the triplet contribution 
to the pn —>■ pnr] cross section in terms of that for pn —>■ drj independent of the 



details of the meson exchanges responsible. Generalising the result of Ref. 
if we parameterise the low energy deuteron cross section estimates as 



apn^d^iQ)^a^Q{l + bQ), (6.70) 
then the integration over phase space leads to the prediction that 



a 



\ b ^f^ 1 X 



2^V 2{l + y^TT^)\ 

(6.71) 
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where x = Qjed- 

Deviations from the predictions of Eq. (|6.7CI| ) in our results will mainly arise 
from corrections to the approximation of Eq. ( |6.69|) in the scattering domain and, 
in particular, from D-state effects. 
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7 Results 



Since the value of the 77-nucleon couphng constant is very poorly known, we shall 
start by neglecting i] exchange, to discuss it later. One would expect, on the 
basis of Eqs. (I^ID and (g^ , that with pure p or TT exchange the pn : pp cross 
section ratio would be a factor of 5, which is not far from the observed value of 
6.5 ±1.0 1^. This agreement is, however, destroyed by the final-state interaction, 
which gives a bigger enhancement for the 1=1 state than for the / = (see 
Fig. 2). Pure p or vr exchange would then predict a ratio of 3 or even less. 




5 10 15 20 25 30 35 40 



Q (MeV) 



Figure 5: Variation of the total cross sections for the pp pprj, pn — > prirj, 
and pn — > dr] reactions with excess energy. The data are taken from references 
[0 (square), (open triangle), (star), (open circle), (closed circle), [0 
(closed triangle), and (inverted triangle). The curves are the model predictions 
with standard parameters, evaluated without rj exchange. 



It was suggested several years ago |]TT1 that the cross section ratio could be 
increased by invoking a destructive interference between p and vr exchange in the 
1 = 1 channel, and hence constructive in the / = 0. With the values of the 
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parameters given in Section 4, we find the p-amplitude to be about three times 
stronger than that of the vr, with the uj term being another factor of three smaller. 
These values then lead to the broad overall agreement with the pp pprj and 
pn —y pnrj total cross sections shown in Fig. 5. The pn : pp ratio is a little too 
large but, if the p coupling were increased by a mere 5%, then this would be 
reduced to the experimental factor of around 6.5 0. 

Whereas the absolute values of the cross sections and the pn:pp ratio depend 
very sensitively upon the parameters, the same is not true for the shape of the 
energy dependence of the pp — ^ ppr] total cross section where the near-threshold 
values need to be enhanced. The ?7-nucleon final state interaction, that we have 
neglected in our analysis, could well supply such an enhancement. The incorpo- 
ration of such an effect is important but not straightforward and we have made 
no attempts in this direction. 

The largest cross section shown in Fig. 5 near threshold is that for pn — > dr], 
but its dominance is mainly a reflection of the faster rise of the two-body phase 
space as compared to those of the three-body ppr] and pnr] final states. To see 
this effect in greater detail, we show in Fig. 6 the ratio of the pn — > pnr] to the 
pn dr] total cross section on a linear scale, where the effect of the rapid phase- 
space variation can be clearly seen. The theoretical prediction of Eq. (|6.71|), based 



upon the extrapolation theorem |]T8|, overestimates the ratio and this would be 
made slightly worse if the small 1 = 1 contribution were also included. The 
situation is changed in the full calculation and this is mainly due to the deuteron 
D-state which increases the np —>■ dr] estimate by up to a factor of two. This 
large effect is due to the deuteron wave function occuring linearly in Eq. (|6.68|) . 

Any residual discrepancies in Fig. 6 may be partially experimental in origin 
since the data were obtained using a deuterium target and the values of Q de- 
duced by kinematic fitting [§, ^ The uncertainties in the pnr] final state are 
particularly large and one cannot rule out a systematic shift of a few MeV. In 
order to get better resolution in Q, it would be an advantage to measure the 
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10 20 30 40 

Q (MeV) 

Figure 6: Ratio of the pn — >■ pnrj and pn dr) total cross sections at the same 
excess Q. Data from Ref. 0] are compared to the full calculation (solid 

curve) and the simplified final-state-interaction approach of Eq. ( |6.71|) (broken 
curve) which only includes 1 = contributions. 



spectator proton directly. Such a procedure is now feasible [|4 



The results shown in Fig. 5 have been obtained after neglecting rj exchange. 
Now, because of the large riN elastic amplitude, pure rj exchange with an rjNN 
coupling of G^/An = 1 would lead to a. pp ^ ppr] cross section twice as big 
as that for n exchange. If the other parameters were not altered, a value of 
G^/Air = 0.25 would reduce the pn:pp ratio prediction to below 5. To reproduce 
this experimental charge ratio with GI^/Att = 1, the p-exchange amplitude would 
have to be reduced by a factor of two, and vector meson dominance normally 
leads to estimates that are far more reliable than that. Since the charge ratio 
is only weakly affected by initial distortion or choice of A^A^ wave functions or 
form factor parameters etc., our model indicates that G^/An < 1 and is probably 
much smaller. 

The pp ppr] cross section is a function of cos^ because the two initial 
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protons are identical. From Eq. ( p.l5| ) and ( [4.36| ), it is seen that to first order 
in the (i-wave amphtudes the differential cross section should be linear in cos^ ft^. 
This is illustrated in Fig. 7, where the predictions, without t] exchange, have been 
increased by 25% before comparing them with the results of [|r^ . 

The agreement with experiment is reasonable, except for the largest angle 
point. The large slope parameter quoted in []T0| is due to the weight of this point 
in the free fit to the data and we see no way whereby such a value could be 
achieved within the current approach. It should also be noted that this point 
affects the overall scale, because the angular distribution is normalised to the 
integrated cross section found in previous studies 0. 



c 













1' ■ ' 








1 



0.4 0.6 



Figure 7: Distribution in the cm. angle of the r] in the pp — >■ pprj reaction at 
Q = 37 MeV (closed circles) and Q = 16 MeV (open circles) from [|r^. The 
predictions of our model with standard parameters without r] exchange, shown 
as the corresponding solid and broken lines respectively, have been increased by 
a factor of 1.25. 



Though small, the slope is still much larger than that found for ^yp ^ rjp 
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or TT^p Tjn [^. It is here important to note that the photoproduction data 
are maximal for 0^ 90° where the pion-production results show a minimum. To 
get agreement with the pn/pp total cross section ratio, the p and vr contributions 
must cancel in the pp case, and this enhances the curvature and leads to the 
agreement shown in Fig. 7. For pn prirj/dr] the interference is constructive and 
we would expect the slope to be an order of magnitude smaller. 

There yet, no data on the analysing power Ay in pp pprj, though 

there is an approved proposal at COSY-11 to attempt the first investigation . 
There are also no measurements of Ay in n^p ^ rjn but, if we consider only 
p exchange, we predict that Ay for pp pprj should be of the same size as 



that for 'yp ^ rjp |]3^ at the same value of the rj-N momentum. With the 



parameterisation of Eqs. (|3.29| , |3.3CI|) , one expects that the analysing power to 
vary as 

Ay = 2 A^*^ sin Or, cos . (7.72) 

The predicted energy dependence of the maximum value of the analysing 
power A™^^ is shown in Fig. 8. The values are small because of the node in the 
input of Eq. ( p.30| ) at r] ^ 0.4, and so they may be sensitive to the neglected 
p-wave and pion-exchange terms. It may be of interest to note, however, that the 
maximum occurs for Q ^ 10 MeV, which is a region where the acceptance of the 
COSY-11 spectrometer is significant ES]. 
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Figure 8: Predicted value of the pp — pprj maximum analysing power coefficient 
of Eq. ( |7. 72| ) as a function of excess energy, estimated assuming only vector-meson 
exchange. 

8 Conclusions 

We have been able to describe, within a one-meson-exchange picture, the main 
features of the experimental data on rj production in nucleon-nucleon collisions 
near threshold. Although most of the energy dependence is reproduced, this is 
mainly an effect of phase space, combined with the nucleon-nucleon final-state 
interaction. The remaining low-Q deviations in both the A^A^ —>■ NNrj and 
pn drj reactions could be due to the ?7-nucleon final-state interaction that was 
not included in the calculation. 

The good agreement with the magnitude of the pp pprj total cross section 
may, in part, be fortuitous given the uncertainties in the the coupling constants 
and other parameters, including those arising from the use of vector-meson dom- 
inance to extract the input amplitudes. Although we have not adjusted the input 
in order to fit the data, it should be stressed that there remains significant flexi- 
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bility in the numerical predictions arising from these uncertainties. 

A more serious worry is the sensitivity to the short-range variation of the 
nucleon-nucleon wave functions arising from the form of the production operators. 
The Paris NN potential , that we have used in all our estimates, is essentially 
local and this has a significant repulsive core that suppresses the wave functions 
at short distances. On the other hand, the Bonn potential ^ has a non-local 
component that allows the core to be much softer and this results in much less 
small-r suppression. This means that the evaluation of the pn — > drj cross section 
on the basis of the integrals in Eqs. (|6.671 , |6.68|) is about three times bigger for 
the Bonn deuteron wave function than that of Paris. The difference would be 
even larger without form factors to regularise the production operators at short 
distances. We would argue that, because we are using local production operators, 
then it is more logical to use wave functions corresponding to a local potential. 
The np —>■ drj cross section is sensitive to the deuteron D-state wave function 
at short distances and the ratio of this to the S-state also differs significantly 
between Paris and Bonn. 

In our approach, the p-exchange amplitude is bigger than that of the vr and it 
is particularly reassuring that the destructive interference required to give good 
agreement for the ratio of the pn pnrj/pp pni] cross sections is such as to 
give the correct shape of angular dependence in the pp pprq case. It should 
be noted that if the vr-exchange term were to dominate, then the slope of the 
differential would have the sign opposite to that of experiment. This illustrates 
the importance of analysing all the //-production channels simultaneously. 

The uncertainty in the rjNN coupling constant makes it impossible to esti- 
mate the //-exchange term with any reliability. It cannot be too large without 
destroying the good agreement with the pn/pp ratio and this allows us to de- 
rive the upper limit within the model of G^^^/47r < 1. It should be noted in 
this context that the rj contribution in semi-phenomenological nucleon-nucleon 
potentials p^, |4l| is to be associated with the exchange quantum numbers and 
not necessarily purely with the true ?7-meson. 
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Although our predictions for the slope in pp pprj are compatible with the 
data, it must be admitted that the experimental measurements are still far from 
definitive in both statistical and systematic aspects. A significant improvement 
is to be expected in the results through the use of the WASA detector . This 
is designed for the study of rare decays of the rj but the production of the meson 
will also be investigated as a by-product under very good conditions. One test of 
our model is the prediction that the rj angular distributions in the proton-neutron 
cases should be much flatter than for proton-proton. This could be checked in the 
pn di] reaction, where the measurement of all fast final-state particles means 
that the kinematics may be well defined even without detecting the spectator 



proton The only measurements extant showed an i] distribution that was 
consistent with isotropy 0. 

Another partial check on our model could come from a measurement of the 
beam analysing power in pp pprj. This is predicted to have a maximum, 
albeit rather small, at quite low values of Q. The predictions are not conclusive 
because of the lack of n^p ^ rjn analysing power input to complement the 
extensive unpolarised differential cross section data which should appear in the 
near future p6[] . 

The gravest drawback in the calculation is the neglect of t] rescattering. There 
are as yet no well-behaved ?7-nucleon potentials in the literature and the short 
range of the r^-production operator means that this effect will depend upon the 
details of the 77-nucleon interaction and not just the on-shell t-matrix. A consis- 
tent three-body treatment of the r] NN final state in rj production is a challenge 
still to be met. 
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